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Abstract —In this paper, we investigate the artificial noise- 
aided jamming design for a transmitter equipped with large 
antenna array in Rician fading channels. We figure out that 
when the number of transmit antennas tends to infinity, whether 
the secrecy outage happens in a Rician channel depends on 
the geometric locations of eavesdroppers. In this light, we first 
define and analytically describe the secrecy outage region (SOR), 
indicating all possible locations of an eavesdropper that can cause 
secrecy outage. After that, the secrecy outage probability (SOP) is 
derived, and a jamming-beneficial range, i.e., the distance range 
of eavesdroppers which enables uniform jamming to reduce the 
SOP, is determined. Then, the optimal power allocation between 
messages and artificial noise is investigated for different scenarios. 
Furthermore, to use the jamming power more efficiently and 
further reduce the SOP, we propose directional jamming that 
generates jamming signals at selected beams (mapped to physical 
angles) only, and power allocation algorithms are proposed for 
the cases with and without the information of the suspicious area, 
i.e., possible locations of eavesdroppers. We further extend the 
discussions to multiuser and multi-cell scenarios. At last, numer¬ 
ical results validate our conclusions and show the effectiveness 
of our proposed jamming power allocation schemes. 

Index Terms —Jamming, massive MIMO, outage, security. 

I. Introduction 

Massive multiple-input multiple-output (MIMO) systems, 
where an enormous number of antennas are deployed at the 
base station, have become a hot research area in recent years 
Q, 0. As the number of antennas goes to infinity, the 
effect of uncorrelated interferences and noises can tend to zero 
asymptotically by using only simple linear transmit/receive 
techniques 0, leading to intensive growth in spectrum and 
power efficiency 0. When used for beamforming, massive 
MIMO leads to sharp beam patterns as well as low power 

Manuscript received December 19, 2014, revised April 21, 2015. The editor 
coordinating the review of this paper and approving it for publication was Z. 
Han. 

Jue Wang, Jemin Lee, and Tony Q. S. Quek are with Singa¬ 
pore University of Technology and Design, Singapore 487372. Email: 
{jue_wang, jemin_lee, tonyquek}@sutd.edu.sg 

Jue Wang is also with School of Electronic and Information Engineering, 
Nantong University, Nantong 226019, China. 

Fanggang Wang is with State Key Laboratory of Rail Traffic Control 
and Safety, Beijing Jiaotong University, Beijing 100044, China. Email: 
wangfg@bjtu.edu.cn 

The contact author is J. Lee. 

This research was partly supported by the Temasek Research Fellowship, 
the National Natural Science Foundation 61401240, 61201201, U1334202, 
the Key Grant Project of Chinese Ministry of Education (No. 313006), the 
Fundamental Research Funds for the Central Universities (2015JBM112), and 
the State Key Lab of Rail Traffic Control and Safety (No. RCS2014ZT09). 

Part of this work has been presented in IEEE Global Commun. Conf., 
Workshop on Massive MIMO: From theory to practice, Austin, TX, Dec. 
2014. 


leakage to unintended directions Q. Due to these attractive 
properties, massive MIMO becomes a promising technique 
for future communication systems such as the fifth generation 
cellular system Jh). In the meanwhile, it can be anticipated that 
massive MIMO will also become crucial in security related 
applications. 

Secure communication in wiretap channels has been stud¬ 
ied for decades since the seminal work 0. Corresponding 
studies have been further extended to different type of wiretap 
channels iU, ID, fading channels ifTOl . ifTTl . MIMO channels 
nm-insi, and networks ESI, Eli. The research topics span 
a wide range from information-theoretical contributions such 
as secrecy capacity analysis and rate region characterization 
to practical transmission design issues including precoding, 
user scheduling, and artificial noise (AN)-aided jamming. For 
a complete review of the most lately approaches, see El, 
ED- Regarding the communication secrecy, the emergence of 
the massive MIMO technique brings new opportunities and 
challenges. Recently, physical layer security techniques using 
massive MIMO have drawn increasing attentions in the litera¬ 
ture. In ll^ - ll2l . the secrecy rate in massive MIMO systems 
has been analyzed using large system analysis and secure 
precoding schemes were designed. In ll24l . ll2l . it has been 
shown that massive MIMO can benefit the detection of active 
eavesdropper who performs attacking on the channel training 
phase ll26ll . Note that the above-mentioned approaches require 
either the channel state information (CSI) of eavesdroppers 
can be known, or their existence can be detected. For the 
scenarios that eavesdroppers are completely passive and their 
CSI is unknown, AN-aided jamming ll27l can be a feasible 
solution. Only recently, the AN-aided jamming approach has 
been applied for massive MIMO systems in ll2l . 1291 and was 
shown to be beneficial for communication secrecy. 

In this paper, we study the secure communication in massive 
MIMO systems via AN-aided jamming. Differently from l2l . 
we consider the scenario that eavesdroppers are randomly 
located around a legitimate transmitter equipped with large 
antenna array, and all channels follow Rician distribution. In 
this case, the geometric locations (described by both the angle 
of arrival and distance to the transmitter) of the legitimate 
receivers and eavesdroppers become essential in the secrecy 
outage analysis, which highlights the main difference between 
our work and l2^ . The motivation of our paper is based on 
the following considerations: 1) Since the beam towards the 
legitimate receiver becomes sharper and the power leakage to 
other directions becomes trivial in massive MIMO systems, 
it is doubtful whether jamming is still beneficial for secrecy. 
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and 2) as the number of antennas grows, the dimension of the 
jamming space increases and jamming power needs to spread 
over a large number of directions, which makes conventional 
uniform jamming inefficient with massive MIMO. Regarding 
these issues, two questions are raised: 

1) Does conventional uniform jamming still benefit the 
secure communication in massive MIMO systems when 
Nt goes to infinity? 

2) Is there more efficient scheme rather than uniform 
jamming in the massive MIMO setup? 

In this paper, we will answer these two questions by making 
the following contributions: 

• For the massive MIMO Rician fading channels, we 
analytically describe the secrecy outage region (SOR) 
as geometric locations of eavesdroppers that can induce 
secrecy outage. The concept of SOR further has been used 
to characterize the secrecy outage probability (SOP). 

> With the information of the suspicious area where eaves¬ 
droppers are possibly located, we derive analytical ex¬ 
pression of the SOP in the presence of one legitimate 
receiver and multiple passive eavesdroppers. After that, 
it is proved that conventional uniform jamming is still 
useful in terms of reducing the SOP when any eaves¬ 
droppers are located within a certain distance range to 
Alice, which we call it as the jamming-beneficial range. 
This conclusion provides an answer to the first question. 

• For uniform jamming, the optimal signal and jamming 
power allocation is investigated for different scenarios. 
We further devise practical directional jamming algo¬ 
rithms, either with or without the information of the sus¬ 
picious area. The proposed directional jamming schemes 
use the jamming power more efficiently to further and 
substantially reduce the SOP, which provides answers to 
the second raised question. 

The rest of this paper is organized as follows: Section II 
provides system model. Section III describes the SOR, further 
provides an analytical expression of SOP and a jamming- 
beneficial range. Optimal jamming power allocation is stud¬ 
ied for uniform jamming in Section IV, and in Section V, 
directional jamming algorithms are proposed. In Section VI, 
the SOR is discussed for multiuser and multi-cell scenarios. 
Section VII concludes this paper. 

II. System Model 

In this section, we first present the network model. As an 
important concept in subsequent analysis, we further define the 
normalized crosstalk between two wireless links and introduce 
its characteristics. Then, the AN-aided secure transmission and 
the definition of SOP are described. 

A. Network Model 

We consider the network shown in Fig. [T] where a transmit¬ 
ter (Alice) equipped with Nt antennas transmits to a single¬ 
antenna user (Bob) in the existence of L external passive 
single-antenna eavesdroppers (Eves 1,..., T). Alice uses beam¬ 
forming for the data transmission to Bob, while jamming 
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Fig. 1: Description of the network layout. 


with AN in other spaces (or directions). We define the set 
of receivers = {b,ei, ...,6^} where b denotes Bob and ej 
{I = 1,...,L) denotes Eve 1. Considering Rician fading, the 
channel between Alice and receiver i is given by 


h, = 


K, 


1 + K, 


1 


l + K, 


gj, Vi e Ir 


( 1 ) 


where Ki is the Rician AT-factor, gi 

fast fading part whose elements follow CN(Q, 1) distribution 
(complex normal distribution with zero mean and unit vari¬ 
ance). Eor uniform linear array with inter-antenna spacing do 
(in wavelength), the line of sight (LOS) component hi can be 
written as the steering vector at incident angle Oi. 


h, = s(6»,) = 

( 2 ) 

where Oi is the LOS angle of receiver i. In addition, we 
consider large scale fading where di is the distance from 
Alice to receiver i, and a is the path loss coefficient. 

We consider a practical scenario that Eves are uniformly 
distributed within an angular range Ae = [dmim ^max] and a 
distance range Vg = [D„iin (0e),D 

max (^e)]5 where Djjji[j(de) 
and Dmax(^*e) are functions of 6e € Ac, defining two borders 
of this area. Throughout this paper, we use 


Ksns = m, dc) I 9e G Ae, 4 G De} 0) 


to define the suspicious area. In practice, if Alice has only lim¬ 
ited information of Dniin(0e) and D„iax(^e). she can assume 
the two boundaries are defined by constant values, d^^in and 
dmax- For instance, if Alice knows nothing about TZsus, she can 
set Ac = [0, 27r], X>e = [0, Tmax], indicating that the suspicious 
area (from Alice’s point of view) spans the entire space with 
radius rmax- The effectiveness of this assumption, referred to 
as “constant boundaries” and defined below, depends on that 
how accurately it can describe the real TZsus- 

Definition 1 (Constant Boundaries): To facilitate practical 
design, it is convenient to set the two boundaries of TZsus 
to be constants such that Dinin(fi*e) = rfmln, Dmax(fi*e) = 
draax, V^e G .Ae. 
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fl-FA(cPo(^)), PVi<4-<l 

h-P ’a (cPo(-^)) + E^(Fa(cP„.i(4-)) -Fa(cP™.2(-^))), PVm+i<7^<PVm 


B. Normalized Crosstalk 


For hi, hj (j € X^) defined as O, the following asymptotic 
results hold as iVt —oo EO): 

1 




— j ^ i. 


(4) 

(5) 


where = denotes the 


approximation that is asymptot- 

A KiKj _J j_ A 

^ {1+Ki){l+Kj) 


ically accurately] Ki-j = “ 

Stemming from (|5]l, we intro¬ 
duce the following definition. 

Definition 2 (Normalized Crosstalk): Define the normal¬ 
ized crosstalk between nodes i,j as 

2 


Si-,j{di,0j) — 


—h^h, 

Nt * ^ 




( 6 ) 


Lemma 1: The normalized crosstalk Si-j(9i,9j) has the 
following characteristics: 

1) Si-j{9i, 9j) is a sinc-like function composed of one main 
lobe and multiple side lobes. 

2) With fixed 9j and random 9i ^ W(0niin, ^max) which is 
uniformly distributed between 9jnin and 0max, the CDF 
of Si-j{9i, 9j) can be written as (|7]i (see top of this page) 
where the definitions of CPq, CPm.i, CPm ,2 and Xa(') 
are referred to Appendix lAl 

3) With fixed 9j, the feasible range of Si-j{9i,9j) is 0 < 

Si-j < where < Ki-j is determined by the 

distribution range of 9i, i.e., A = [6*min, ^'max]- 

Proof: See Appendix lAl ■ 


allocation coefficient as 


p. 

A jam 

Ptnt 


Z^n=l 

Pu 


(9) 


For ease of description, we assume that Bob and all Eves 
share the same noise covariance being Nq. Moreover, we 
consider maximum ratio transmission (MRT) for precoding of 
the data symbol Xh, i.e., w;, = In this case, according 

to ([8]l, the SINR at receiver i is given by 


SINK,; = 


Pb^nhfwbi 


^0+d-“Ell^n|hfv, 


( 10 ) 


We assume that the Eves are not colluding, but consider the 
most-capable Eve, which has the maximum receive SINR, to 
define the secrecy rate as 


Rl = [log 2 (l + SINRb) - log 2 (l + SINRe.„,ax)]+ (H) 

where SINRe,max — maxSINRe, and [x\^ = max{a;,0}. We 
say a secrecy outage occurs if is less than a target rate 
Pth, hence the SOP is defined as 


Pout = Pr{P^ < Pth}. (12) 


III. Secrecy Outage Analysis 

In this section, we first introduce the secrecy outage region 
(SOR) which describes all possible locations of Eves who can 
cause secrecy outage. Analytical expression of the SOR is 
derived for uniform jamming, then the SOP is studied with 
variant shapes of TZsus- At last, a jamming-beneficial range 
is derived to show that uniform jamming is still useful in 
reducing the SOP. 


C. Secrecy Transmission Scheme 

We use linear precoding for data transmission, while AN 
symbols s„ are sent in the space defined by Vn,n = 1, 
to degrade the channels of Eves. Eor the null space-based 
jamming lIZTl . it holds that N = Nt — I and v„ S null(hb). 
The received signal at receiver i is given by 

Vi = \J Pbrf~°hfWbCCb + y^ v„s„-fnj, Vi £ P 

n—1 

( 8 ) 

where Wb £ is the precoder for Bob, a:b is the unit- 

norm data symbol, and rii is the additive Gaussian noise. 
Moreover, Pb and P„ respectively are the powers allocated 
to Bob and the n-th^amming direction, with total power 
constraint such as X]n=i = Ptot — Ph where Ptot is the 
total available transmit power. We define the jamming power 

* In this paper, we focus on the large antenna regime, and will use equalities 
instead of approximations for brevity. 


A. Secrecy Outage Region 

In the large antenna regime, all fast fading effects are 
completely averaged out as shown in (@1) and Q. Therefore, 
whether the secrecy outage occurs or not, will be essentially 
determined by the geometric location of Eve. In this light, we 
introduce the SOR defined in the following. 

Definition 3 (Secrecy Outage Region): The SOR is defined 
in terms of polar coordinates as 

7^SOR = |(0e,de) I 1™ P^ < Pth 1 ■ (13) 

At-too J 

Herein, we note that P^ is a function of 9e and de- 

In the large antenna regime, secrecy outage occurs if there 
exists at least one Eve within the SOR. If all Eves locate 
outside of the SOR, the target secrecy rate Pth can be 
guaranteed. To characterize the SOR, we first evaluate the 
received SINRs assuming uniform jamming. 
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Csi^) 


C2{^) 

Ci{4>) 


(l + (1 - - 2«‘») Aot<^ 

(l + (1 - ct>)Ptotd^‘^Nt - Aot<^ + (1 - ■ 


( 22 ) 


Lemma 2: With uniform jamming in null(hb), the SINRs 
at Bob and any Eve can be, respectively, written a^ 

SINR^j = (14) 

SiNRuj ^ - Arfe“lVtSe;b(6>e) - 

l + de-“Pjam(l-Se;b(0e)) 

where, and hereafter, we use the notations A = ^, ^jam = 
for brevity. Note that in ( fTST i, Se;b(^*e) is the normalized 
crosstalk between Eve and Bob as defined in (|6]). Considering 
fixed 0b, we hereafter write Se;b as a function of only 0^. 

Proof: Since span(v„) = null(hb), jamming causes no 
interference at Bob. Applying (IHi to (fTOl i. we get (fT4l i. Noting 
that A^ = — 1 and A = (flOl) . we have 


SINR^j = 


‘Ihfwbl^ 


lVo + de-“fel:Ve^W„|2 


n—1 


(16) 


By applying (|6l), the numerator of (fTbl) can be written 
as P\^df°‘NtSe-h{de)- On the other hand, noting that Wb 
and v„,n = l,...,Wt — 1 constitute a complete orthog¬ 
onal basis of the A^t-dimensional vector space, we have 

|heW„p= ||he|p - |hoWbp = Wj(1-Se;b(^^e)) in the 

denominator of (fTbl) . Therefore, (fTSl l can be obtained. ■ 
Remark 1: The result in (fl4l i leads to a constraint on f 
(defined in (|9])), written as 


(j) < = 1 - 


2flth _ 1 

Ptotdf'^Nt 


(17) 


which stems from the fact that the jamming power cannot be 
too large, otherwise, even without Eves, the target rate i?th 
cannot be guaranteed since the remained signaling power is 
too small. Unless otherwise specified, we assume (EJ can 
always hold via proper power allocation. 

Erom Lemma |2] we characterize the SOR for the uniform 
jamming as follows. 

Proposition 1: With uniform jamming in null(hb) and 
given (j), the SOR is described as 


^SOr('/>) = S (^e, de) I de < 0e), Se;b (de) > 


(18) 


(19) 


where 

d'"l(</),de) = (C'i((^)Se;b (de) " C2(</>)) “ 

Clif) = - (1 - 4')PtotNt2^^'' - p 

l + (l-</.)Ptotdr7V,-2«- 

C2{f)=Ptotf (21) 

and C' 3 ((/)) is given by (l22l i shown at the top of this page. 


^Hereafter, if one notation is applied for any Eve, we will use the subscript 
e instead of ei for the sake of brevity. 


Proof: Substituting (fT4l i and (flST l into (fTTT i and using the 
definition of SOR in (fTSl l. we can obtain the value of dgl {f, dg) 
in (O. Note that the value of dgpfjOe) should be positive, 
this straightforwardly introduces the constraint on the mini¬ 
mum value of Se;b(de) such as Se;b (de) > C^{4)) where C‘i{4i) 
can be readily obtained by letting C'i((/))se;b(de) — C 2 {(t)) > 0. 


Erom Lemma[T] Se;b(de) in ( fT9] l is a function with one main 
lobe and multiple side lobes, resulting in a multi-lobe shaped 
SOR. In order to gain some insights from this complex shape 
(as will be shown in the simulations), we focus on several 
critical security-related metrics as 

« Largest radius of the main lobe ( dJ^'^Q) and the m-th side 
lobe (d"j^).- de;'o “d d^j„ can be obtained by replacing 
Se;b(de) in (fl^ respectively with ATeib and ATeibPVm 
(defined as (l54l i in the proof of Lemma [TJ, such that 

d2o^{CPm.-,h-C2m- (23) 

d'"^^ = (Ci(<^)iTe;bPV™ - C2(</>)) “ . (24) 

Since dg-^Q is much larger than f 0, dg-^Q can be 

considered as the largest distance of the SOR. Eor any 
Eve whose distance to Alice is larger than dg'^Q, we can 
conclude that it causes no secrecy outage regardless of 
its LOS direction. 

« Largest angle difference ASmax of the SOR: Eor any 
Eve whose angle difference to the LOS direction of Bob 
is larger than we can conclude that it causes 

no secrecy outage regardless of its distance to Alice. If 
Se;b(de) < Ca((/>), V0e > de, we Can Write 


A0 


max 



(25) 


Clearly, to have smaller SOR, we expect to reduce both dg'^Q 
and A0„iax- To minimize d^-’g in ( l2Tl l. we need to minimize 
Ci{(j)) while maximizing C' 2 (</>); however, since C^{4)) = 
, this results in a maximized (correspondingly, 

a larger A0max), which is not desired. It is clear that a trade¬ 
off in (p exists in balancing the effects of both dg^g and 
ASiuaxj which can be formulated as jamming power allocation 
problems, as described in the following sections. 

At last, by setting (/) = 0 in Proposition [T] we obtain the 
following corollary: 

Corollary 1: Without jamming, i.e., 0 = 0, the SOR can 
be found as 


^SOR — 


(de, dg) I de < 


/ AotiVt2^^»-5e;b (de) \ 

Vl + Ptotd,-“iVt-2^*i> ) 



where the superscript (•)"•* stands for “no jamming”. 

Proof: The corollary is directly obtained by setting 0 = 0 
in Proposition [T] ■ 
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Differently from (fTST l, the constraint on Se;b(^e) vanishes in 
( |26] |, indicating that the SOR now is extended to the entire 
angular domain. Moreover, compared with ( l2?t . we see that 
de,o in no-jamming case0on the contrary, is reduced compared 
to uniform jamming. In conclusion, uniform jamming induces 
two opposite effects: the beneficial one is that the SOR can 
be squeezed in angular domain, and the disadvantage is that 
the SOR is enlarged in Bob’s direction, i.e., the main lobe. 
Illustration of the SOR changing caused by jamming will be 
shown later in simulations. 


B. SOP Analysis 

With a single Eve uniformly distributed in 'R-sns, using the 
derived SOR, the SOP is given by 


outjSingleEve 


Area(7esoR(<^) n7^ sus) 

Area (7^sus) 


(27) 


where Area( ) denotes the area of a certain geometric region. 
Considering that there are L Eves uniformly distributed in 
Ttsns, the SOP of the entire network can be written as 


T^out — 1 (1 7^out,singleEve) 


(28) 


Prom dZTl) . the SOP is determined by the overlapping area 
between two geometrical regions. If 7 ?.sor(<(')07^sus = 0, zero 
SOP is achieved. Recalling Q, as well as (l2?t and (l25l l. two 
sufficient conditions of 'R-'50'r{4>) C T^sus = 0 can be written 
as 


(fT2li as Vout = 1 — Pr{i?^ > i?th} and recall (fTTTi . we have 


Pr{7?b > 7?th} — P'siNR";'^^,, 


/ 1 -f SINR^j 

1 



(31) 


where Agjj^j^uj (•) is the CDP of which is 

given by = (^sinrJ^ (a^)) since all Eves are 

independently distributed. Using (fTSl l. we have 


'^SINRj'’(^) — Pi' S 'Se;b(^e) ^ 


Pi. 


Pv.Nt 




(32) 


jam , 


where both 6^. and de are random. Since 9^ and de are 
independent, (l32l i can be presented as 


PsiNRe(a:) = 


P. 


- jam 




+ Pa 


fdSz)<iz. (33) 


where fde{z) = — is the PDP of corresponding to 

max min 

the uniform distribution between two boundaries defined by 
Ve = [dmin, dmax]- Then, Pout IS directly obtained as (O. ■ 
Practically, (l30l l can be used for jamming power alloca¬ 
tion. As stated in Remark 1, Alice can arbitrarily adjust the 
value of the constant boundaries in the design, based on the 
information about Pgus that she has. Particularly, if Alice 
knows nothing about Psus (i-C-, she assumes Ae = [0,27r] 
and Pe = [0,rmax]), minimizing Pout becomes equivalent to 
minimizing Area (PsOR (<(>))■ 


or D„,ax(0e) = 0, V |0e - 0b| < A^max (29) 

where de,o is defined in (|2^ and A^max is in (|25] |. The 
physical insight of ( |29] | is clear: when an Eve is far away 
or its angle difference to 0b is large, it does not cause outage. 

Por the general case with arbitrary shape of TZsus, Pout 
in (l27l i can be numerically evaluated and further applied to 
jamming power allocation design. However, due to the non¬ 
regular shapes of Psor( 0) and TZsus, closed-form expressions 
of Area (Psor(<^) n Psus) as well as the SOP in (l28T l do 
not exist for the general case. Yet, by considering constant 
boundaries of TZsus as described in Definition 1, (l27l i can be 
written in an integral form as the following proposition. 

Proposition 2: With constant boundaries of TZsus, i e., 
Ac = [0min,0max] and Pe = [dmin, dmax], and Uniform 
jamming in null(hb), the SOP can be given as 


Pout — 1 ~ 


Fs. 


1 p 

^ -I 1£ 


PbAft2«th 


P. 




2z 


d^ - (F ■ 

^max '^mm 


dz \ . (30) 


where Pie b(') defined in (I7]i. 

Proof: Por the ease of analytical description, herein we 
utilize the CDP of the normalized crosstalk in (l7]i. Pirst, rewrite 


^Hereafter, for the conditions where the couesponding notation can be 
applied for either the uniform-jamming or no-jamming cases, we ignore the 
superscript “uj” or “nj” for brevity. 


C. Jamming-beneficial Range 

Based on Proposition |2] we find a jamming-beneficial range 
defined in dmax (i e-, the larger constant distance boundary of 
TZsus) as follows. 

Proposition 3: A constraint on dmax that makes the uniform 
jamming beneficial in reducing the SOP is given by 


< 


PtotiVt2««' 


^e;b 


1 + Ptotdr“iVt - 2««> 


(34) 


where is the largest feasible crosstalk value defined in 
Pemma J] 

Proof: See Appendix IbI ■ 

Remark 2: Proposition [3 shows that when Eves are located 
close enough to Alice, uniform jamming is always beneficial 
in reducing the SOP. Clearly, this range expands with larger 
dh, as well as larger or larger Pth- On the other hand, 
the range shrinks with larger Nt or Ptot- 

Moreover, we note that 01 has a similar form of that 
described for the SOR without jamming, i.e., PgoR in (l26l l. 
Recalling the definitions of the largest distance of SOR in (l23t 
and (l2l . the physical insight of Proposition[3]can be explained 
as follows: as long as TZsus Cl Pg’oR f 0, there always 
exists an optimal with which the SOP can be reduced by 
uniform jamming, compared with the SOP without jamming. 
The optimization of is discussed in the next section. 


IW. Jamming Power Allocation 

In this section, considering uniform jamming, we investigate 
the optimal jamming power allocation that minimizes the SOP. 
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The problem can be simply described as 

min Pout, s.t. 0 < 0 < 1. (35) 

<l> 

In practice, Alice may have different accuracy levels of 
information about TZsus as follows: 

1) Alice knows nothing about the suspicious area, or only 
partial information about the suspicious area such as Ae 
only (or only); and 

2) Alice knows exact information about the suspicious area, 
i.e., both Ae and Ve- 

For these two cases, we respectively investigate the jamming 
power allocation in the following. 


A. Jamming with None/Partial Information about TZsus 

When Alice knows nothing about TZsus, minimizing the SOP 
becomes equivalent to minimizing the area of SOR, which can 
be calculated as 

p‘ 2 'k -I 

Area(p"V(<^)) =y^ - 0e))'d0e (36) 

where 0e) is dehned in (fT^ . Note that P-sor(()>) is 

composed of many side lobes. We use to denote 

the m-th side lobe, and PsORi('^) denote a group of 
side lobes with indices described by the set X. For the case 
that Alice knows Ae or Ve, we can simplify the problem by 
minimizing partial, other than the entire area of 'Pgo]^(()>) such 
as 

Area (PsV.i'('/')) = ('^soR,mW) (37) 

m^X' 

where X' is the set of the concerned side lobe indices, 
determined by either Ae or T>e- Using (l3^ (or (ITTT i) along 
with (fT^ . the areas can be numerically calculated and the 
optimal (j) can be easily founded via one dimensional linear 
search. Since it is difficult to derive closed-form expression 
for Area evaluate the area of PgoR,m('/') 

in the following corollary for a special case to further provide 
some discussions. 

Corollary 2 (Area o/PgQj^ ^(i^)).' With = 0 and the 
path loss coefficient being a = 2 (which corresponds the 
free space propagation Area ^PgoR ^ can be upper 

bounded as 

Area(K^'oR,„(«) < jjA - 2a(i.)) 

(38) 

where Ci{(f) and C 2 {(t>) are dehned in (l20l i and (1211 1. 

Proof: See Appendix O ■ 

Remark 3: In (l38] l. it is shown that the area of every side 
lobe is inversely proportional to Nt, indicating that the SOR 
side lobes can asymptotically vanish with ultimately large Nt- 
Moreover, it is inversely proportional to nA, which means 
that the area of the SOR will rapidly decrease for the side 
lobes with large indices, i.e., with large angle difference to 
0b. This result indicates that Eves from different directions 
(i.e., within different side lobes) have different signihcance in 


causing secrecy outage, hence should be treated differently in 
the jamming design. 


B. Jamming with Exact Information of TZsus 

With the information of TZsus, Alice can calculate and apply 
the value of 'Pout in the design (at least numerically)!! using 
either (|28] | or (l30l l. Although in practice, (l35l) can be readily 
solved by one dimensional linear search, it fails to provide 
the optimal f in closed form. In the following corollary, we 
provide closed-form solutions and discussions for a special 
case. 

Corollary 3 (Jamming power allocation for given 9e): For 
constant boundaries of TZsus and given 0e, which is equal for 
all Eves, the optimal f can be determined as 


where 


= 


min{(/)7\</'o} ^oe[0,l] 


(39) 



<^o 


df'^NtPtot 

Ptot j _ 

l-EPtotdb°Aft-2-^th min 

(1 — Se-,h{de))Ptot 


(40) 

(41) 


Proof: See Appendix iDl ■ 

Note when the optimal jamming power 

decreases with dg and Se;b(0e), whereas it will increase when 
(^°P^ = (pQ. The part that dominates the hnal result in (|39] | 
depends on the value of 0e- Detailed discussions will be 
provided in Section VII along with simulations. 


V. Directional Jamming 

In this section, we propose directional jamming algorithms 
to allocate jamming power more efficiently than uniform 
jamming, based on the following facts: 

1) With the information of TZsus, Alice can perform jam¬ 
ming only to the suspicious directions instead of the 
entire null space of hg. 

2) Without information of TZsus, jamming towards different 
directions also needs to be treated differently, as stated 
in Remark [3] 

At a cost of slightly increasing the implementation complexity 
compared with uniform jamming, directional jamming is able 
to substantially reduce the SOR In following subsections, we 
present power allocation algorithms for directional jamming 
with and without the information of T^sus- 


A. Directional Jamming with the Information of TZsus 

When jamming is not uniformly performed, from (flOl i. the 
SINR at Eve is represented as 

JK“|hfwbp 

No + de-“hf Vdiag (p) V^he 

'^The calculation requires the knowledge of T)e and Ae- Clearly, uniform 
jamming is not optimal in this condition. However, for the ease of analysis, 
we first devise the optimal power allocation for uniform jamming; then, the 
resulted jamming power can be allocated directionally to further improve 
efficiency. 


SINR;^j 
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Algorithm 1 Directional jamming with the information of 

"^sus 

1: Initialization: Update the information of 0b, and T^sus 
at Alice. 

2: Assuming null space-based uniform jamming, hnd 

= argmin 'Pout 
4> 

through one dimensional linear search over cj) G [0,1). 

3: Select a iV-dimensional subspace iVn,n = and 

N < Nt) from nulljhb} according to (l46l l. where 0v„ is 
dehned as dH). 

4: Equally allocate = (jP'^^Ptot to the selected beams 

popt 

in step 3 such as 


where V G i^NtxN matrix that spans the jamming space, 
with the n-th column vector being v„, and p = (Pi,P/v) , 
where P„ is the power allocated to the n-th jamming direction 
as dehned in (|8]l. Correspondingly, the SOR now can be 
described as 


^sor(p) 


{(de,0e) |de <de'(P,0e)} 

Y 2^^>>(l-0)PtotiVtSe;b(0e) 

_Vl + (l-<^)Ptotd,-“7V,_2«th 

- s(0e)^Vdiag(p)V"s(0e)) “ 


(43) 


(44) 


where p = §(•) was dehned in (|2]). The superscript 

stands for “directional jamming”. 

Design directional jamming using (l42l i induces high com¬ 
plexity especially when N is large, since changing any ele¬ 
ment in the iV-dimensional vector p requires re-calculation 
of P-sor(p). Hence, we alternatively propose a two-step 
suboptimal power allocation method for directional jamming 
in Algorithm [T] which hrstly hnd the optimal jamming power 
assuming uniform jamming, then reallocate it directionally 
based on a criterion of jamming subspace selection. In Step 
2, Pout can be calculated numerically using (|28] | or (l30l l. 
depending on the available information of Psus- Note that 
0b, dh and Psus are long term parameters, thus the updating 
period of Step 1 can be much longer than the computation 
time required by the other steps. 

Since Pgus is dehned by physical angles, it is necessary to 
set up a mapping between the jamming space and physical 
angle to concentrate the jamming power towards Psus- We 
propose a heuristic subspace selection method for Step 3. First, 
map v„ to physical angle as 

0v„ = arg max |s(0)^v„P. (45) 

^^[”"2 ’"2 ] 

After that, Pj°^ is equally reallocated to the beams whose 
indices are 

Af = {n I 0v„ G A}- (46) 

popt 

The power allocated to each beam is now In practice, 

v„ is not necessarily in null(hb) and an alternative is to hnd 


Algorithm 2 Iterative directional jamming without informa¬ 
tion of Psus 

1: Initialization: Update the information of 0b and db at 
Alice; set an initial jamming power allocation vector po, 
satisfying ||po|| < Ptot'P’^^^', set the initial iteration index 

j = 1 - 

2 : for n = 1 to TV do 
3: Update 

Pj{n)= argmin Area 7^s■'oR( [Pj(1)> 

xG[0,iCinax) y V 

Pj{n - + 1), ...,Pj-i{N)] 

n— 1 N 

where PA^) “ L Pj-AA- 

2=1 2=n+l 

4: end for 

5: if IIpj - Pj-ill > £ then 
J = J + 1; go to step 2. 

7: else if IIpj — pj_ i||< e then 
8 : return 

9: end if 



v„ as the column vectors of a iVt-dimensional DFT matrix for 
the following reasons; 1) selected columns of the DFT matrix 
can form a good substitute of null(hb), as Nt ^ oo 1331; 
2) using pre-dehned DFT basis as the jamming space avoids 
channel inverse calculation, which induces high computation 
complexity especially when Nt is large 12^ : and 3) most 
importantly, the structure of DFT matrix provides very sharp 
beam pattern towards the physical angle 0v„ in ( |45] |. therefore, 
the beam selection criterion (l46l l can be very efficient since 
with sharper beams, there will be less jamming power leaked 
outside of TZsus- 


B. Directional Jamming without Information of 'R-sus 

Without any information of 'R-sus, the objective of direc¬ 
tional jamming power allocation becomes to minimize the area 
of T^goR(p) in (l43t for 0e G [0, 27r], which is calculated as 

P‘ 2 'k -| 

Aiea(n%^{p)j = -(dfJ(p,0e)) d0e. (47) 

A general closed-form expression of WT\ is not available, 
and its convexity is unknown. Hence, numerically minimizing 
Area (n'^AR(p)) is NP-hard. To overcome this, we propose 
Algorithm 2, which iteratively hnds the optimal n-th element 
of p while keeping the others hxed. 

Algorithm |2] provides a sub-optimal solution which reduces 
the complexity by degrading the original problem to one¬ 
dimensional linear search. However, for large Nt, the complex¬ 
ity is still huge since during each main iteration, N ^ 0{Nt) 
times of linear searching are required to fully update p. Hence, 
Algorithm|2]is not suitable for some scenarios where 0b and dh 
change fast. In this light, we propose a simplihed algorithm. 
Algorithm 3, to further reduce the complexity. In Step 2 of 
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Algorithm 3 Simplified directional jamming without informa¬ 
tion of TZsus 

1: Initialization: Update the information of 0b and db at 
Alice; Initialize Areamin = Inf-i 
2: Calculate Se;b(^m), 171=1,M. Determine m[,m 2 , sat¬ 
isfying that Se;b(0m;) > Se-h{dm'J > Se;b(0m),VTO ^ m[ 
and m ^ 

3: for ^ = 0 to 1 do 

4: Find 

Pboundary = Urg min Area (pboundary)^ • 

5: Calculate Areag^o^ (pboundary) according to gT]); 

6: if Areag^oR (pb^^ndary) < Areamin then 

7: Area„,i„ = Areag^o^ (pboundary) ; 

9: else 

10 : continue 

11 : end if 

12 : end for 


Algorithm 3, 9m is the mean angle of the m-th side lobe of 
^sor(p) (117 = 0 denotes the main lobe). In Step 4, Pboundary 
follows the structure such as 


Pboundary — { 0; 5 0; •■• / ; 


(48) 


where m'l and are the indices of the two dominating side 
lobes, which are located most closely to the main lobe (from 
both sides). The derivation of Algorithm 3 is described in 
Appendix |E] 

With Algorithm 3, jamming is performed in only two dom¬ 
inating directions in the neighborhood of 0b, for the reasons 
that 1) for the region with large angle difference to Bob, 
allocating much jamming power is inefficient since 72.gQ]^ in 
this region is generally very small; 2) for the directions highly 
in-line with 0b, jamming should be avoided as it will cause 
severe interference to Bob. Note that for every realization of 
(p, only single time of linear search is required in Step 4. 
The complexity is irrelative to N (which is large in general), 
hence can be greatly reduced compared to Algorithm |2] As 
a possible extension, more than two dominating directions 
can be involved in the design while the trade-off between 
complexity and performance exists. 


VI. Extension to Multiuser and Multi-cell 
Scenarios 


A. Multiuser Transmission 

When multiple legitimate users (i.e., multiple Bobs) are 
presented in massive MIMO systems for Rician channels, 
the multiuser interference between Bobs is trivial as long 
as their LOS angles have large difference, which can be 
readily ensured via user scheduling. On the contrary, the 
multiuser interference to Eves can be seen as equivalent 
jamming considering that single-user decoder is adopted at 
Eves, which is likely to happen when Eves are low-cost 
devices. Hence, when multiuser beamforming is applied for 
Bobs, the received multiuser interference at Eve becomes equal 
to the directional jamming, transmitted towards other Bobs’ 
directions. Consequently, the SOR of an objective Bob will 
be shrunk in the directions of the other Bobs. Denote the set 
of all legitimate users as = {bi, ...,bf/}. Similar to (l43T l 
and (l44l i. which describe the SOR for directional jamming, the 
SOR of user b^ S in the presence of multiple users can 
now be described as 

^“o'R'’“(PbJ = {{de, 0e) I de < , de)} (49) 




(Pb„,de) 


/ 2^^Hl-</>)AotA^tSe;b„(ge) 


- s(0e)^Vdiag(pb„)V^s(0e)^ 


(50) 


where V = [Wb, Vj] S spans the equivalent jamming 

space, in which Wb = [wb,^,] ,it' = ^ u 

spans the signaling space for the other legitimate users, 
while Vj G null (Wb) is the jamming space. Corre¬ 
spondingly, the power allocation vector can be divided into 


two parts as pb„ 


~ T 

pL 


"T 

uPb„ 


where Pb„,sig = 


T 


-Pbi, 


,Pl 


bu 


is the signal power alloca¬ 


tion vector for all legitimate users except for bu. Eor 


fixed pb„.sig, in order to minimize 


Area ^ 


7? 

/\'< 


MU.b„^- 
SOR 



user hu, the directional jamming algorithms, i.e.. Algorithm 2 
and 3, can be directly applied herein. 

Considering communication secrecy for the entire multiuser 
transmission system, the optimization problem can be reason¬ 
ably re-formulated as a min-max problem such as 


min max Area (Pb j) , 

Pbu bu \ / (51) 

s.t. ||pb„||i < Ptot,hu 

The main challenge in solving ( fSTT i is that allocating power for 
one user affects the SORs of other users. Hence, the power 
needs to be jointly allocated for all users, and the complexity 
of such joint optimization can be very high. Hence, it is 
desirable to develop simplified algorithms for the multiuser 
scenario. 


We focus on the single-cell and single-user scenario in 
previous sections. In this section, we now show how the 
SOR can be affected by multiple users and cells. We also 
provide discussions on the design of secure transmission in 
these scenarios with future research challenges. 


B. Multi-cell Network 

In multi-cell massive MIMO networks, it is commonly 
assumed that the training pilots are reused among cells. 
Correspondingly, pilot contamination results in imperfect CSI 
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Fig. 2; Radius of the main lobe/side lobes of TJ-sor vs. (j). 
Parameters setting: 0b = 0°, db = 100m, Nt = 100, and 
i?th = lObps/Hz. 

estimation as well as nonnegligible multi-cell interference 
from Alices in adjacent cells. Denote 9^ as the angular 
direction of Bob in Cell i seen from Alice in Cell j, and Cell 
0 as the objective cell where Bob 0 exists. Some major effects 
of imperfect CSI and multi-cell interference on the SOR of 
Bob 0 are described as follows: 

• Due to pilot contamination from Bobs i, \/i ^ 0, the SOR 
of Bob 0 will be enlarged in the directions of 0°^^. 

• Multi-cell interference to Bob 0 will isotropically enlarge 
his SOR. On the other hand. Eves in Cell 0 are equiva¬ 
lently jammed by the multi-cell interference from Alice j 
in cell j, Vj ^ 0, especially in the directions of 0p and 
0j jj. Thereby, the SOR in these directions can be shrunk 
and the shape can be non-continuous in these regions. 

A general analytical description of the SOR in the multi¬ 
cell network is challenging, since it is determined not only by 
the network topology, but also by the locations of all pilot¬ 
contaminating users in adjacent cells. Moreover, the compli¬ 
cated shape of the SOR makes difficulties in calculating and 
minimizing the corresponding area. Nevertheless, in practice, 
pilot scheduling and reuse schemes can be utilized to alleviate 
these adverse effects which are caused by pilot contamination, 
e.g., m, GS). 

VII. Simulation Results 

Simulation results are shown in this section. We set do = 
0.5, a = 3, Ptot = IW, Nq = 10“®mW and for simplicity, 
assume strong LOS environment such that Ke-\^ —>■ 1. In 
this parameter setting, the receive SNR is 20dB when the 
transmitter-receiver distance is 100m. 

At first, using (|2^ and (l24l i. Fig. |2]provides a description of 
^sor('^) terms of the radius of the main lobe/side lobes. 
It is shown that the radius of the main lobe is monotonously 
increasing with (j), indicating that reducing the signal power to¬ 
wards 0b enlarges the SOP in this direction. Clearly, allocating 
additional jamming power to the direction of 0b will further 
enlarge this radius, which suggests that jamming directly in the 



Fig. 3: SOP vs. power allocation coefficient <f). Parameters 
setting: dh = 100m, 0b = 0°, De = [50m, 100m], = 

[-15°, 15°], L = 10, and i?th = lObps/Hz. 


legitimate user’s direction should be avoided. This conclusion 
coincides with the concept we followed for the design of 
Algorithm 3. 

Moreover, as (p increases, the radius of the second side lobe 
is reduced first and then increases after a certain value, e.g., 
(p sa 0.7, indicating an optimal jamming power allocation in 
terms of minimizing the SOP in this direction. The side lobes 
with index m > 3 can be completely eliminated with proper 
jamming. The results show that we can design jamming based 
on the partial information of TZsus- For example, in Fig. |2] if 
we know that Eves are located in the direction ranges of side 
lobes with indices larger than 3, then, allocating p = 0.2 is 
enough to secure the communication and the remaining power 
can be allocated to data transmission. 

Next, we depict the SOP vs. p in Fig. [3 with randomly 
generated 0e and d^ within the range = [—15°,15°] and 
De = [50m, 100m], respectively. Note that the smoothless of 
the curves is not due to the lack of simulation trials, but caused 
by the fact that Pout is a piecewise function, as shown in (j?]) 
and (l30l) . In addition, the SOP rapidly increases to 1 when 
p exceeds p’^^^ in (Eli. From Fig. [3 we can see that 1) 
the SOP with optimal p is much smaller than that without 
jamming, i.e., (/) = 0, as anticipated in Proposition [3 2) the 
SOP decreases as Nt increases since larger Nt results in higher 
received power at Bob and less leakage to Eve. Moreover, the 
optimal p increases with Nt because with larger Nt, allocating 
more power to data transmission is not efficient in increasing 
the achievable rate of Bob because of the logarithmic slope 
of the rate function; and 3) the SOP is further substantially 
reduced with directional jamming. 

Fig.g] shows the optimal jamming power coefficient p°'^'^ 
as a function of the normalized crosstalk Se;b, for uniform 
jamming. We compare the derived p°'^'^ in ( |39] | with Monte 
Carlo simulations and show a good match between them. 
From Fig. @1 we first observe that each curve is divided into 
two parts, respectively representing that po or pg^'^ dominates 
the optimal result in ( |39] |. The division is emphasized using 
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Fig. 4: Optimal jamming power coefficient vs. Se;b under 
Corollary 3 with X>e = [50m, 350m] and iith = 5 bps/Hz. 



Fig. 5; Area of the secrecy outage region vs. c?b, results are 
shown for variant jamming algorithms. Parameters setting; 
Nt = 50, Rtb = 5, and 6»b = 0°. 


a vertical dash line for the case with dh = 150m. When 
db = 100m, in the -dominating region, the curves with 
Nt = 50 and 100m coincide with each other since Nt does not 
affect 00 in (l4Tli . In the -dominating region, 0°?^ increases 
with Nt with given Se;b- The region-division in Fig. |4] can be 
explained as follows; in the left region, Se;b is small enough 
so the channels between Eves and Bob can be considered as 
asymptotically orthogonal, i.e., ho G null(hb). In this case, 
more signal power is leaked to Eve with larger Se;b, hence we 
need more jamming power allocated in null(hb) to degrade 
Eve’s channel. However, in the 0g^^-dominating region, the 
value of Se;b is large, which indicates that channels from Alice 
to Eves and Bob could be highly aligned, i.e., he ^ null(hb). 
In this case, jamming in null(hb) is not efficient to degrade 
Eve’s channel and the jamming can be a waste of transmit 
power. Thus, the optimal jamming power starts decreasing 
with Se;b- 

In Eig. |5] for the case that without information of TZsus, we 
compare the area of SOR achieved by different algorithms. Eor 



Eig. 6; Secrecy outage probability vs. d^, results are shown 
for no jamming and uniform/directional jamming algorithms. 
Parameters setting; Nf = 100, i?th = lObps/Hz, 0b = 0°, 
L = 10, Ae = [-30°,30°], and Ve = [50m, 200m]. 


uniform jamming, the optimal 0 is found via one dimensional 
linear search by minimizing the area described in ( |36] |. Eor all 
curves, the SOR enlarges with increasing dg. since larger dg 
results in weaker signal power received at Bob. With uniform 
jamming, the area of the SOR can be reduced approximately 
by half compared with the non-jamming case. This consider¬ 
able reduction of SOR area together with the least implemen¬ 
tation complexity make uniform jamming still a good option 
for practical system design. Using the directional jamming 
with Algorithm 2, the SOR area can be further reduced, but it 
induces the highest complexity among all schemes. At last, we 
note that for the directional jamming. Algorithm 3 achieves 
slightly larger area of SOR than Algorithm 2, where the 
difference becomes smaller especially for small db. However, 
the implementation complexity can be greatly reduced by 
Algorithm 3, which makes it being a reasonable choice that 
strikes a compromise between complexity and performance. 
It can also be confirmed from Eig. |5] that even without any 
information of TZsus, directional jamming can still be utilized 
to enhance communication secrecy. 

We depict the SOP vs. db in Pig. |6] for a particular example 
scenario, where 7?.sus is defined by Ac = [—30°,30°] and 
Ve = [50m, 200m]. Pour schemes are compared, i.e., without 
jamming, uniform jamming, and directional jamming schemes 
Algorithm 1 and Algorithm 3. Correspondingly, we also depict 
the SORs achieved by each schemes in Pig.|7]to help us better 
understanding the relation between the SOR and SOP. We note 
that in Pig. |6] the smoothless of the curves is caused by the 
area calculation of the intersection between two complicated¬ 
shaped regions, not by the lacking of simulation trials. 

As shown in Pig. |6] all jamming schemes can achieve 
lower SOP compared to that without jamming. In the small-dg 
region, directional jamming schemes outperform uniform jam¬ 
ming. However, the performance improvement from uniform 
jamming to directional jamming Algorithm 1 is small because 
Algorithm 1 directionally allocates jamming power towards 
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(a) Without jamming 



(b) Uniform jamming 
200 


(c) Directional jamming (Algo. 1) 


90 200 



(d) Directional jamming (Algo. 3) 

90 200 



Fig. 7: Description of the corresponding SORs for the 4 
schemes in Fig. |6] when db = 150m. 


TZsus to suppress the SOR side lobes (as shown in Fig. |7] 
(c)), whereas when Bob is close to Alice, most side lobes 
are already very small with uniform jamming. Moreover, in 
the small-db region, directional jamming Algorithm 3 achieves 
the lowest SOP for the reason that it focuses on only two 
dominating SOR side lobes (as shown in Fig. |7] (d)) hence the 
jamming power can be used more efficiently in Algorithm 3 
than Algorithm 1, where jamming is uniformly performed to 
all directions within TZsus- In this example, the two dominating 
side lobes are covered by TZsus, thereby they contribute more 
in the SOP calculation compared with the other side lobes. 
However, in a different scenario where TZsus does not cover 
the two dominating side lobes, it cannot be concluded that 
Algorithm 3 always outperforms Algorithm 1. 

At last, we note that as dh increases, the performance 
of Algorithm 3 degrades and Algorithm 1 outperforms the 
others. The reason is, when dh is large, all SOR side lobes 
correspondingly become large as shown in Fig. |7] (d). In this 
case, besides the two dominating side lobes, the impact from 
the other side lobes cannot be simply ignored as that has been 
done in Algorithm 3. In conclusion, in practice, appropriate 
jamming scheme should be determined according to both the 
available information of TZsus and the location information of 
Bob such as dh- 


VIII. Conclusions 

In this paper, we first define and analytically describe the 
SOR for secure communication in massive MIMO Rician 
channels, and derive expressions of the SOP. We then de¬ 
termine a jamming-beneficial range, indicating that uniform 
jamming is useful in reducing the SOP when the distance 
from Eve to Alice is less than a threshold. Optimal jam¬ 
ming and signal power allocation is investigated for uniform 
jamming, furthermore, for both conditions with and without 
the information of the suspicious area, we propose directional 
jamming algorithms, which makes use of the jamming power 
more efficiently and further reduce the SOP. We further extend 



Fig. 8; Description of s(x), the cross points and peak values. 


the discussions to multiuser and multi-cell scenarios where 
future challenges are also described. In conclusion, we claim 
that uniform jamming still helps the communication secrecy 
in massive MIMO systems, and the proposed directional 
jamming outperforms conventional uniform jamming schemes. 


Appendix A 
Proof of Lemma[T] 

Proof of T): We can rewrite Si-j{di,dj) in (|6]l as a function 
of Ai-j as follows 


9j) 


Ki-jS {Ai-j). 


(52) 


By applying fST] (14)] to U-j in (|6]l, s{x) in (l52l i can be 
represented as 


s(a::) = 


1 , 

1 sin^(7V£7r(ia;) 
N'^ sin^(7rdx) ’ 


x = 0 

X ^ 0 


(53) 


which is a sinc-like function, has one main lobe and side lobes 
with decreasing amplitudes. 

Proof of 2): In order to describe the CDF of Si-j, we first 
characterize s(a;) in (l5^ by some cross points and peak values, 
shown in Fig. and defined in the following. 

Definition 4: The cross points and peak values are defined 
as 


« Peak Values: The peak value of the m-th side lobe of 
s{x) can be approximated by 


PV 


m 



1 

7r2 {m + \f 


(54) 


which is obtained by noting that sin(a;) « x when x is 
small, and it becomes asymptotically exact as Nt is large. 
PVq = 1 corresponds to the main lobe. 

• Cross Points: When u < PVm, CPm,i{u),i = 1,2 
denotes the i-th cross point between y = u and y = s{x) 
in side lobe m (m < M), CPo(m) is the cross point in 
the main lobe. 

Using (l52] |. we have Fs-.^{Ki-ju) = Pr{s(Ai;j) < m}. 

According to Fig. [8] Pr{s (Ai-j) < u} can be evaluated by 
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calculating the probability that falls within the discrete 
intervals determined by GPm,i{u) and CPo(u). On the other 
hand, recalling that A^y = |sin0i — sin0j| and Oi follows 
uniform distribution, the CDF of A^y can be written as 




min (sin ^ (min (1, 2 ; + sin0j)) , 


— max (sin ^ (max (—1, —z + sin0j)), . (55) 


Using Fa{z) to describe the probability that A^y falls within 
the described intervals leads to (ITJ. 

Proof of 3): Given Oj, the feasible range of A^y is deter¬ 
mined by Ai, i.e., the angle range that node i distributed in. 
Use to denote this feasible range (shown in Fig. |8), it 
holds that C [0,1]. Given Ai, the feasible range of s^y 
can be determined as [Q,K^jj3\], where f3\. = max s{x). 

It is clear that < 1, leading us to the conclusion. 


Appendix B 

Proof of Proposition^ 


For the ease of description, we first define 


h (x) = 


1 -f xd7°^Nt - 2«ti* ■ 


(56) 


Letting t = 


z +Pjam _ hlPA, and to = — where 

<^1+-^ jam 


0-2 = h{Ptot), we rewrite 


as 


'Pout — 1 ~ 


^•5e;b (^) ^ 


2z 


^2 - (P 

^max ^min 


dz } . (57) 


Since F'se b(’) ™ increasing function, T^out decreases with 

t in the domain of Therefore, if the following condi¬ 

tions are satisfied, we can conclude that jamming is beneficial. 

1) There exists a positive value of Pjam, which holds t — 
to > 0 for any z S 2?e; and 

2) z GVe such that to < 

where condition 1) defines the scenario, in which jamming 
can always result in a larger t than to (corresponding to the 
no-jamming case), thereby leading to lower T^out since it is 
a decreasing function in t. Moreover, condition 2) ensures 
that to is less than the maximum feasible value of Se;b, 
otherwise, T^out is always zero either with or without jamming 
under condition 1), hence the benefits of jamming cannot be 
concluded. For 1), it is equivalent to prove that 


-- > - 

Ol + Pjara 02 


(58) 


has a positive solution of Pjam- Note that with proper param¬ 
eter setting that f is no larger than described in ( fTTI i. 

we have ai > 0 and 02 > 0, and note that if 


d 


a 

max 


< 02 


(59) 


then z“ < 02 holds for any z. In this case, (|58j can be 
equivalently rewritten as 


Tjam ^ 


(oi - 02 ) 2 :° 
02 — 2 :“ 


(60) 


Noting that < 0, it is clear from (l60l l that Pjam > 0, 

thus ( |59] | is a sufficient condition to satisfy 1). On the other 
hand, as to < a sufficient condition that satisfies 2) is 


d 


a 

max 


^ „max 

< Se;b “2- 


(61) 


Here, according to Lemma [T] we have < 1. Hence, the 
intersection of (l5^ and (IhTl i is equal to (ImT i. and recalling the 
definition of 02 leads to the final result. 


Appendix C 

Proof of Corollary[2] 


First, from (ES, m, Se-o{0e) in the m-th side lobe can be 
presented as 


Se;b(^e) 


Ke-h sin^ (Ntndx) 
sin^ {irdx) 


m m + 1 
7^’ Ntd 


(62) 


In (l6^ . sine functions appear in both the numerator and 
denominator, making it infeasible to obtain closed-form results 
in further analysis. Hence, we find an upper bound of (l62l i by 
fixing the value of x to the left end point of its domain, i.e., 
X = -ffh, in the denominator. 

I\tci 


Se;b(^'e) < sin^ {Ntirdx). 


(63) 


For the ease of description, we consider only the condition 
that Oe > 0. With Oo = 0, we can rewrite x (defined as a; = 
I sin^e — sin 0bI) in terms of Og as 


X = sin 9g < 9 g. 


(64) 


Note that the side lobes that are close to the main lobe are more 
important in contributing to the area of the SOR, as a result, the 
value of 9g that should be concerned is very small. Therefore, 
the upper bound in dMl) can be very tight. Combining (l6^ 
and ( l64l l we have 

Se;b(0e) < sin^ {Ntnd9e) . (65) 

Substitute dbSl) into (fl^ . and calculate the area by integral, 
we get 


Area 


1. 


< 


1 

G,4m 2 
1 

eGAm 2 


{Cyf)sg.oi9g)-C2md9g 


( 66 ) 


sin^ {Ntirddg) - C2{4>) ) d0e 


(67) 




Nfird 


( 0 ) “ C '2 (<())) dp 


1 


dNfd V Tr'^rrP 


,Cyf)-2C2{f) 


( 68 ) 

(69) 
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where Am = the physical angle range of the 

TO-th side lobe. To obtain ( |66] |. we use a = 2, and is 
obtained using ( |65] |. From ( |67] ) to dhSl l. we use the variable 
substitution rj = NtirdOg, then ( |69] | is obtained by simple 
integral calculation of elementary functions. 


Appendix D 

Proof of Corollary[3] 

When 0ei = 0e, V/ = Tout in (l30l l is determined 

only by such that 

\ L 


'Pont — 1 ~ 
where 


2z 


/ 

/max[(io(0,^e),<^niin] ^max ^min 


dz 


(70) 


do((/), 9e) = ((1 - <^)Aot) + 4'Ptot'^ Se-,h{9e) “ 4>Ptot 

= [gm- (71) 


is the distance threshold. Given 6^, any Eve with a distance 
to Alice smaller than do{(j),6e) can cause secrecy outage. 
Noting this, dTOll calculates the overall outage probability 
by assuming that de follows uniform distribution. Clearly, if 
mmdo{(j),9e) < dmin, the integral range in (iTOl) becomes 

[dniinj dniax] hence a zero SOP can be achieved; otherwise 

if min do 9g) > dmax, a definite outage occurs with proba- 

(j) 

bility Rewrite g{(j)) in (17X1 1 as 

= Se.\y{9g)gi{(j)) + (1 - Se;b(de))fl2((^) (72) 

where gi{(()) = d((l - (fyPtot), g 2 {(t>) = -4>Ptot and h{-) 
is defined in (l56l l. Clearly, minimizing (iTOl i is equivalent to 
minimizing (Gall. Since pi((^) is concave and g 2 (</>) is linear, 
g{(j)) is concave. Hence, letting = 0, the optimal cj) 

satisfying cj) < (as in (flTl O is given by in dXOl) . 

Note that if do(0g^^,de) < dmm, setting is 

not the only choice since the solution of do((/), de) = dmin (if 
exists), denoted as (po, also achieves zero SOP. Therefore, we 
need to check the value of po and compare it with (p°g^^ to 
determine the final optimal jamming power allocation. Note 
that in (l72l i. the value of gi{(j)) is usually much less than 
g 2 {'P)- Hence, when Se;b(de) is not so large, g{(j)) can be well 
approximated by a linear function as g((/)) ss Se;b(de)d(iXot) — 
(1 — Se-h{9e))Ptot'P- Using this approximation, we get po = 

ge;b(^e)^(Ftot)~C^niin 

(1 —Se;b(®e))Aot 

If po ^ [0,1], it is not a feasible solution in practice. If 
po & [0,1], both po and pg^^ are able to achieve zero SOP. In 
this case, we choose the smaller one between pl^'^ and po to 
save jamming power such that p°'^'^ = minj^g*’'^, ^o}- 


Appendix E 

Derivation of Algorithm 3 

Algorithm 3 is proposed based on the following assumptions 
and approximations; 

^To facilitate the analysis, we assume the definite outage will not happen 

by letting dmax to be large enough such that dmax > maxmindo(<^, (^e). 

Se </> 


1) Assuming that the angle ranges occupied by every side 
lobe (and the main lobe) of 72.g^Qf(^(p) are approximately 
the same (which is assuming that there are in total 
one main lobe and M side lobes within the angle range 
[^, ^]), an upper bound of the integral in (l47l i can be 
approximated by 


2 ) 


M 

AreauB (7^s'oR(p)) ^ {df{p,9m))^ 


(73) 

where the area of every side lobe is upper bounded by 
the area of its enclosing sector. 

We assume that s{9m)^'Vm = 1 and s(d„)^Vm = 
0,Vn 7 ^ m. In practice, by defining v^, = 
this assumption is easy to realize with large iVj, where 
asymptotic orthogonality holds. Now, according to (l44l) . 
dtKP,0m) in (l73l) can be written as 


di\v,9m) = 


(^^m Prn^ 


(74) 


and Pm = ^ 


where n - 
Where a,„ - 

is the m-th element of p, with Pm being the jamming 
power allocated on the m-th side lobe. 


Erom (l73l l and (l74l i. and given fixed p, we rewrite the original 
area minimization problem as 


M 

min ^ (^Qm 

m—1 



2 


M 

S.t. Pm = pPtot,0 < Pm < a.. 

m—0 


(75) 

(76) 


By checking the Hessian matrix, it is easy to show that the 
objective function in (TTSI) is concave. To minimize a concave 
function, clearly, the optimal solution can be found only 
on the boundaries of the domain defined by (|7^ . Recalling 
that the area of the side lobes decreases rapidly with larger 
lobe index, and being aware that jamming should be avoided 
within the main lobe to prevent degrading Bob’s channel, we 
further simplify the problem by checking the boundary of the 
domain as described in dXHl l. According to these discussions. 
Algorithm 3 is obtained. 
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